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Abstract 



Assuming that p -mesons exist in a quark-gluon plasma at temperatures 
close to the QCD phase transition, we calculate the dilepton production rate 
from q-q annihilation via a p -meson state using Vector Meson Dominance. 
The result is compared to the rates from direct q-q annihilation and from 
7r + -7r~ annihilation. Furthermore we discuss the suppression of low mass 
dileptons if the quarks assume an effective mass in the quark-gluon plasma. 
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Dileptons are one of the most promising signatures for medium effects in the fireball in 
relativistic heavy ion collisions. At SPS there is an indication for an enhancement of the 
dilepton production at invariant masses between 200 and 800 MeV compared to all known 
sources neglecting medium effects 0. Possible explanations of this enhancement are the 
broadening of the width of the p -meson by scattering in the medium and the reduction of 
the p -mass by the on-set of the restoration of chiral symmetry || . 

Furthermore dileptons may indicate the formation of a quark-gluon plasma (QGP) phase, 
in particular at RHIC and LHC 0. In a perturbative calculation the lowest order contribu- 
tion to the dilepton production from the QGP comes from the direct q-q annihilation into 
a virtual photon (Born term) At high invariant masses M this contribution is believed 
to dominate. For low invariant masses of the order of gT, on the other hand, a^-corrections 
become important . Close to the critical temperature g might be as large as 6 || . Hence 
a s -corrections could be important even at invariant masses up to M = 1 GeV or larger. 

In order to avoid infrared singularities and gauge dependent results, HTL resummed 
propagators and vertices || have to be used in perturbative calculations of the dilepton 
production from the QGP [§,0. The dilepton rate follows from the imaginary part of the 
finite temperature photon self energy. Unfortunately, it turned out that 2- and probably 
even 3-loop contributions |TD[ to the photon self energy are of the same order in a s as 
the 1-loop contribution and exceed the latter one even by factors of 3 and more. Thus the 
application of finite temperature perturbation theory suffers from the facts that higher order 
loop corrections are important even in the weak coupling limit and moreover that realistic 
values of the strong coupling constant are not small. 

The importance of non-perturbative effects in a QGP at temperatures within reach of 
heavy ion collisions has been realized in QCD lattice calculations, where for example the 
equation of state, condensates, and hadronic correlators below and above the phase transition 



have been studied [ |TT|j . 



It has been noticed that the equation of state found on the lattice can be described 
perfectly by an ideal gas of massive quarks and gluons (quasiparticles), where an effective, 
temperature dependent mass of the order gT for the partons has been introduced ||. The 
effective quark and gluon mass at temperatures around the critical one T c is of the order 
0.5 GeV. This will lead to a complete suppression of dileptons from the QGP below about 
M = 1 GeV as we will discuss at the end of the present paper. 

The existence of a gluon condensate above T c can be used to construct an effective quark 
propagator in the QGP |T2|. The quark dispersion relation following from this propagator 



has two branches, of which one shows a minimum at finite momentum. This leads to an 
interesting structure of the dilepton spectrum, exhibiting peaks (Van Hove singularities) and 
gaps H[[3|. 



In the present paper we want to study the consequences of the existence of hadronic 
states in the QGP at temperatures around T c as indicated by lattice calculations [|TJ| and the 
gauged linear sigma model [15) on the dilepton production from the QGP. This temperature 
regime plays an important role in relativistic heavy ion collisons as the mixed phase with 
T = T c may exist for a long period in the fireball and might contribute significantly to the 
photon and dilepton production (see e.g. [^6[ ). In particular, we consider a possible co- 
existence of quarks and p -mesonic states in the QGP, because p -mesons are an important 
source for dileptons. Although one expects that p -mesons vanish above T c quickly with 
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increasing temperature fli"4jl , they might be present at temperatures around T c |15[ . 



We assume a simple coupling of the quarks to p-mesonic states in the QGP. Besides by 
direct q-q annihilation electron-positron pairs can then be produced by the process qq — > p — > 
e + e~ , where the decay of the p will be described using Vector Meson Dominance (VMD). The 
result will be compared to the rates following from direct annihilation of quarks (Born term) 
and from ti + tt~ annihilation using VMD. Furthermore we assume a vanishing baryon density 
(zero quark chemical potential), equilibrium distributions for the quarks and p -mesons, and 
the bare mass m p = 770 MeV for the p -meson. 

For the interaction of quarks with p -mesons we adopt the following Lagrangian: 



C>=-- A PlvfC+~rnlp a lx ff a + q [i 7/1 0" - m q + g ^ f < q , (1) 



where p a = d p p^ — d u p p with the p -meson field p p and the quark field q. Here a denotes 
the isospin or flavor index and r a the corresponding isospin matrix. For the quark mass m q 
we will investigate two cases, namely a vanishing bare mass and an effective quark mass 
m g = 0.5 GeV. 

To get an idea about the size of the coupling constant g describing the strength of 
the quark-p coupling we proceed in the following way: By integrating out the p -mesons 
from Lagrangian (fl|) one obtains in lowest order of the derivative expansion the four-quark 
interaction term 



A-quark = ~~ (qiJ^q) , ( 2 ) 



which might be compared with the corresponding interaction term from the extended NJL 
model -G2 (qifiTaq) 2 Hl7| , [l8| . This suggests the identification g = J &m?G2 which serves to 



determine g by taking G 2 from the literature p~T|JT8|] : g m 5 — 6. In the following we choose 



g = 6. Of course, it is not clear whether our choice for mass and coupling of a vector meson 
excitation in a hot QGP is justified. However, in lack of any first principle calculation we 
decided to choose values as suggested by measurements or calculations for the vacuum case. 

The dilepton production rate (here for massless e + -e _ pairs) can be calculated from the 
imaginary part of the photon self energy according to [19| 

dN a 1 ImU^(P) 



d 4 xd 4 p 12vr 4 e E / T - 1 M 2 ' 

where a = e 2 /4ir. Here we use the notation P = (E,p) and p = \p |. 

Using VMD the photon self energy is related to the p° -meson propagator by 

ImIi^P) = - 2 m i p ImD^P) 1 (4) 
9 P 



where g p = 6.07 ||19|| . The most general ansatz for the trace of the imaginary part of the full 



p -meson propagator reads 

„. , ImF 2 ImG 

^ } ~ (M 2 -m? p - ReF) 2 + (ImF) 2 (M 2 - m 2 - ReG) 2 + (ImG) 2 ' y } 



3 



where F(p ,p) and G(po,p) are the longitudinal and transverse parts of the p°-self energy. 
To lowest order the p°-self energy is calculated from the one-loop diagram containing a quark 
loop. It is given by the one-loop photon self energy containing an electron loop multiplied 
by an factor 3/2, where the factor 3 comes from the number of quark colors in the loop and 
the factor 1/2 from the flavor coefficient tr(ro/4), which counts the number of quark flavors 
(u and d). Furthermore the electron charge e has to be replaced by g. 

The one-loop photon self energy at finite temperature can be calculated analytically 
in the high temperature or equivalently hard thermal loop limit However, in this 

approximation there is no imaginary part for timelike photons, M 2 = p 2 , — p 2 > 0, resulting 
in a vanishing dilepton production. Going beyond the hard thermal loop approximation, 
integral expressions for the matter part of the one-loop photon self energy can be derived 
[0 : 



ReF 



ImF 
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ImG 
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(6) 



where uj\ 



k 2 + m 2 n F (uj k ) = l/[exp(u k /T) + 1], and 



1 / 4m 2 

2 l^i^M 2 - 



2 I l'<> w 1 



Am 2 
M 2 



+ p 



(7) 



The second term in the imaginary part, being proportional to 8(—M 2 ), describes Lan- 
dau damping which takes place only for spacelike momenta, i.e. in scattering processes. 
Therefore it does not contribute to the dilepton production. In the hard thermal loop limit, 
Po,p -C k, the first term of ImF and ImG vanishes, whereas the second one reduces to the 
well known Landau damping contribution in this limit 
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We do not take into account the vacuum part of the p self energy because the decay of a 
p-meson into a free quark- ant iquark pair does not take place in the vacuum. The Lagrangian 
(|I|) can be considered as an effective Lagrangian valid only at temperatures at or above the 
phase transition. 

Combining (|]) to (|) we obtain the dilepton production rate (p-quark rate) by numerical 
integration over the loop momentum k (magnitude of the three-momentum). In Fig.l this 
rate is compared to the Born rate (direct quark-antiquark annihilation) [|J and the 7r + -7r~ 
annihilation rate via VMD (p-7r rate) calculated by Gale and Kapusta |19|]. In contrast to 
the Born rate the p-quark rate as well as the p-7r rate show a clear peak in the vicinity of the 
p-mass. For M < 0.4 GeV the p-quark rate agrees well with the Born rate, which both have 
no mass cut assuming vanishing u- and <i-quark masses. For M > 0.4 GeV the p-quark rate 
agrees better with the p-7r rate, which is zero for M < 0.28 GeV due to the finite 7r-mass. 
The agreement of the absolute values of these rates might be accidental as, for example, the 
p-quark rate is enhanced compared to the p-ir rate by the larger number of quark degrees of 
freedom than n degrees of freedom, but reduced by the neglect of vacuum contributions in 
the first one. 

We have chosen a temperature of T = 0.15 GeV in accordance with the critical tem- 
perature for the QCD phase transition as predicted by QCD lattice theory, since we expect 
p -mesons to co-exist with quarks only close to T c . For the momentum of the virtual photon 
we have chosen p = 0.2 GeV, because CERES-SPS data show that the dilepton enhancement 
occurs at small transverse momenta 0.2 GeV< px < 0.5 GeV 0. 

In Fig. 2 the p-quark rate is shown for different temperatures (T = 0.15 GeV and T = 0.2 
GeV). In Fig. 3 the p-dependence of this rate is presented by comparing the rate at p = 0.2 
GeV and p = 1 GeV. 

Finally, we discuss the possibility of an effective, temperature dependent quark mass as 
indicated by comparing the equation of state of an ideal gas of massive quarks and gluons 
with lattice calculations. In the temperature regime between T c and about 2T C , quark masses 
of the order of m = 500 GeV are required to match the quasiparticle equation of state to 
lattice results ||. In Fig.4 the influence of such an effective quark mass on the Born rate is 
shown. The rate vanishes for M < 2m q but approaches quickly the bare mass rate above 
M = 1 GeV. The same behavior can be observed in the p-quark rate (Fig. 5), where the 
p -peak now is absent in the case of finite quark masses. 

If this simple picture for the QGP were true, there would be no dilepton production 
from the QGP below the p-peak, where the SPS dilepton enhancement has been observed. 
On the other hand, although the simple assumption of an effective quark mass might be 
sufficient to explain the equation of state, the correct quark dispersion relation in a QGP 
may look completely different. Perturbative as well as non-perturbative approaches indicate 
two quark branches starting from the same effective mass at zero momentum. The lower 
branch, corresponding to a collective quark mode (plasmino) that is absent in vacuum, 
shows a minimum at finite momentum. The splitting of the two collective quark modes 
and the minimum of the plasmino branch give rise to interesting structures, namely peaks 
(van Hove singularities) and gaps, in the dilepton rate. In particular, low mass dileptons 
are possible coming from electromagnetic transitions from the upper to the lower branch 
P,[13|. Furthermore, a finite width of the quarks from scattering in the QGP will change the 
dilepton production rate in addition. 
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Summarizing, we have calculated the dilepton production rate from quark- ant iquark 
annihilation assuming the co-existence of quarks and p -mesons around the QCD phase 
transition as indicated by lattice QCD and gauged linear sigma model calculations. Using 
VMD the dilepton production rate has been calculated in a similar way as the rate from 
7r + -7r~ rate, replacing the 7r-loop by a quark loop in the latter. For low invariant masses 
the rate turns out to agree well with the Born rate, whereas for higher invariant masses the 
agreement with the p-7r rate is better. Although the coincidence in the absolute value of 
these rates may be accidental, the agreement of the shape is easy to understand (no mass 
cut for massless quarks, p-peak in VMD). 

Concerning SPS data the contribution from the p-quark rate to the dilepton spectrum is 
by far too small to explain the observed enhancement Q . For the p-quark rate is of the same 
order as the Born rate, which leads to a dilepton yield that is about two orders of magnitude 
smaller than the observed one according to hydrodynamical calculations [0. However, this 
situation will be different at RHIC and LHC where a higher initial temperature and a longer 
lifetime of the QGP are expected. 

Finally, we have shown that an effective quark mass as indicated by lattice calculations 
of the equation of state of the QGP leads to a suppression of dileptons with invariant 
masses below about 1 GeV. However, the simple picture of non-interacting massive quarks 
is probably too oversimplified for predicting the dilepton production from the QGP. 



6 



REFERENCES 



[1] G. Agakichiev et al., CERES collaboration, Phys. Rev. Lett. 75 (1995) 1272; N. Masera 
for the HELIOS-3 collaboration, Nucl. Phys. A 590 (1995) 93c; A. Drees for the CERES 
collaboration, Nucl. Phys. A 630 (1998) 449c; G. Agakichiev et al., CERES collabora- 
tion, Phys. Lett. B 422 (1998) 405. 

[2] G.E. Brown and M. Rho, Phys. Rev. Lett. 66 (1991) 2720; B. Friman and H.J. Pirner, 
Nucl. Phys. A 617 (1997) 496; R. Rapp, G. Chanfray, and J. Wambach, Nucl. Phys. 
A 617 (1997) 472; F. Klingl, N. Kaiser, and W. Weise, Nucl. Phys. A 624 (1997) 527; 
W. Peters, M. Post, H. Lenske, S. Leupold, and U. Mosel, Nucl. Phys. A 632 (1998) 
109; W. Cassing, E.L. Bratkovskaya, R. Rapp, and J. Wambach, Phys. Rev. C 57 (1998) 
916. 

PV. Ruuskanen, Nucl. Phys. A 544 (1992) 169c. 

J. Cleymans, J. Fingberg, and K. Redlich, Phys. Rev. D 35 (1987) 2153. 
T. Altherr and P.V. Ruuskanen, Nucl. Phys B 380 (1992) 377; J. Cleymans and I. Dadic, 
Phys. Rev. D 47 (1993), 160; J. Cleymans, I. Dadic, and J. Joubert, Phys. Rev. D 49 
(1994) 230. 

E. Braaten, R.D. Pisarski, and T.C. Yuan, Phys. Rev. Lett. 66 (1991) 2183; S.M.H. 
Wong, Z. Phys. C 53 (1992) 465. 

M.H. Thoma and C.T. Traxler, Phys. Rev. D 56 (1997) 198. 

A. Peshier, B. Kampfer, O.P Pavlenko, and G. Soff, Phys. Lett. B 337 (1994) 235; A. 
Peshier, B. Kampfer, O.P. Pavlenko, and G. Soff, Phys. Rev. D 54 (1996) 2399. 
E. Braaten and R.D. Pisarski, Nucl. Phys. B 337 (1990) 569. 

P. Aurenche, F. Gelis, R. Kobes, and H. Zaraket, Phys. Rev. D 58 (1998) 085003; P. 
Aurenche, F. Gelis, R. Kobes, and H. Zaraket, [hep-ph/9903307 . 
G. Boyd et al., Nucl. Phys. B 469 (1996) 419. 

A. Schafer and M.H. Thoma, [hcp-ph/9811364 , Phys. Lett. B (in press). 
M.G. Mustafa, A. Schafer and M.H. Thoma, in preparation. 
G. Boyd, S. Gupta, F. Karsch, and E. Laermann, Z. Phys. C 64 (1994) 331. 
R.D. Pisarski, Nucl. Phys. A 590 (1995) 553c; R.D. Pisarski, Phys. Rev. D 52 (1996) 
3773. 

J. Sollfrank et al., Phys. Rev. C 55 (1997) 392. 

A. Polled, R.A. Broglia, P.M. Pizzochero, and N.N. Scoccola, Z. Phys. A 357 (1997) 
325. 

Y.B. He, J. Hiifner, S.P Klevansky, and P. Rehberg, Nucl. Phys. A 630 (1998) 719. 
C. Gale and J.I. Kapusta, Nucl. Phys. B 357 (1991) 65. 

V.V. Klimov, Sov. Phys. JETP 55 (1982) 199; H.A. Weldon, Phys. Rev. D 26 (1982) 
1394. 

[21] J.I. Kapusta, Finite Temperature Field Theory (Cambridge University press, New York, 
1989) p.71. 



7 
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FIG. 1. p-quark rate (solid line), Born rate (dashed line), and p-ir rate (dotted line) at T = 0.15 
GeV, p = 0.2 GeV and m q = 0. 
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FIG. 2. p-quark rate at T = 0.15 GeV (solid line), T = 0.2 GeV (dashed line), p = 0.2 GeV, 
and m a = 0. 
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FIG. 3. p-quark rate at T = 0.15 GeV, p = 0.2 GeV (solid line), p = 1 GeV (dashed line), and 



m q = 0. 
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FIG. 4. Born rate at T = 0.15 GeV, p = 0.2 GeV, m q = (solid line), and m q = 0.5 GeV 
(dashed line). 
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